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时 灌 模 糊 奇 异 摄 动 系统 万 滤波 * 


陈 金 香 ， 杨 卫 东 
(北京 科技 大 学 信息 工程 学 院 ， 北 京 100083) 
摘 要: 针对 一 类 变 时 灌 非 线性 奇异 摄 动 系统 ， 本 文 提 出 了 基于 T-S 模糊 模型 的 建 模 方法 ， 并 设计 了 模 
W Ho 滤波 器 ， 通 过 求解 一 组 与 摄 动 参数 < 无 关 的 线性 矩阵 不 等 式 ， 获 得 其 增益 ， 避 免 了 由 s 引 
起 数值 求解 的 病态 问题 ， 所 获得 滤波 器 使 闭环 系统 渐进 稳定 旦 达到 给 定 的 五 w 性 能 指标 。 该 方 
法 适用 于 标准 和 非 标准 非 线性 奇异 摄 动 系统 ， 仿 真实 例 说 明 该 方法 的 有 效 性 。 
关键 词 : AH; NF; Hoo 滤波 ; 线性 矩阵 不 等 式 (LMIs); T-S 模型 


分 类 号 : AMS(2000) 93C40 中 图 分 类 号 : TP13 文献 标识 码 : A 
1 引言 


多 时 标 现象 广泛 存在 于 复杂 电子 电路 、 亲 性 机 嚣 人、 通信 网 络 等 各 种 物理 、 化 学 及 人 工 系 
统 中 。 作 为 处 理由 多 时 标 现象 引发 病态 问题 的 理想 工具 ， 奇 异 摄 动 技术 得 到 了 深入 研究 D7 引 。 
但 目前 传统 方法 研究 非 线性 奇异 摄 动 系统 9， 对 系统 结构 的 假设 还 较 多 。 近 来 ， 模 糊 控 制 理 
论 与 非 线性 奇异 摄 动 系 统 相 结 合 的 研究 ， 较 好 地 利用 线性 奇异 摄 动 系统 相关 理论 ， 解 决 了 非 时 
湾 模 糊 奇 异 摄 动 系 统 分 析 与 控制 问题 {01!9， 但 时 灌 模 糊 奇 异 摄 动 系 统 的 研究 尚 不 多 见 (11 。 文 
HK [16] 利用 LMI 方 法 研究 了 时 滞 模 糊 奇异 摄 动 系统 的 稳定 性 条 件 和 五 状态 反馈 控制 器 问题 。 
本 文 建立 连续 时 滞 模 糊 奇 异 摄 动 系统 模型 ， 并 采用 MI 方法 研究 了 一 类 时 滞 模 糊 奇 异 报 动 系统 
不 依赖 于 摄 动 参数 的 有 Hw 滤波 器 设计 问题 。 该 方法 无 需 快 慢 分 解 且 适用 于 标准 和 非 标准 非 线性 
奇异 摄 动 系统 ， 为 研究 非 线 性 时 滞 奇 异 摄 动 系统 的 互 co 控制 提供 新 思路 。 


2 “时 滞 模 糊 奇异 摄 动 系统 模型 
在 本 节 中 ， 我 们 参考 文献 [16]， 建 立时 滞 模 糊 奇异 摄 动 系统 模型 ， 并 研究 时 灌 非 线性 奇异 摄 
动 系统 的 Hoo 滤波 问题 。 
时 灌 模 糊 奇 异 摄 动 系统 第 i (1 < i 7) 条 规则 为 如 下 形式 。 
Plant Rule i: If & (t) is Fii, and --- and €,(t) is Fy, then 
E(e)a(t) = Aiz(t) + Agiz(t — r(t)) + Buw(t), 
z(t) = Cuz(t), (1) 
y(t) = Cax(t) + Daiw(t), 
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了 全 0 x(t) 
E(e) 一 | , x(t) = | | , 
0 Elmxm £2(t) 


且 zi(t) € R”, a(t) €e RAR, RRAWE, w(t) e R? 为 噪声 信号 (包括 过 程 和 测量 噪 
Æ), y(t) c R* 为 测量 输出 ，z(t) € R! 为 待 估计 的 信号 向 量 ，Fij (7 = 1,2,… ,9) 为 模糊 集 
合 ,， 7 为 规则 数目 ，A;, Adis Bui, Cii, Cai, D21i HEME RERE, Elt) ,g(t) 为 可 测 变 
Eo r(t) < o AJARAN, Hi < 8 < 1， 其 中 7 和 BB 均 为 已 知 常数 。0 < < x1 Wk 
动 参数 。 

采用 标准 的 模糊 推理 方法 一 即 单 点 模糊 化 ， 乘 积 推理 ， 和 加 权 平 均 清 晰 化 ， 得 全 局 模糊 模 
型 为 


其 中 


E(e)az(t) = A(u)x(t) + Aa(u)x(t — 7(t)) + Bi(w)w(t), (0) = 0， 
2(t) = C(u)x(t), (2) 
y(t) = Co(u)x(t) + Dai (u) w(t), 
其 中 
9 
wi = wile), wi(€)) = EO wew) = TT Fu), 
j=1 


wi(é(t) j 
=1 
s(t) ERCH, zalt—7(t)) eR"+™, Alu) = >》 midi, 
i=1 
Tr T r 
Aalu) =Ñ Aa, Bil) =Ñ mBun, Ol) =Y iO, 
i=1 i=1 i=1 


Calu) = uiCn, Dalu) = WiDr E(e) = | k | . 
i=1 i=l 0 elmxm 
对 于 系统 (2)， 如 果 应 用 常规 模糊 系统 理论 和 Lyapunov 方 法 ， 可 得 出 一 组 依赖 于 摄 动 
参数 < 的 稳定 性 条 件 。 但 由 于 s < 1， 系 统 矩 阵 会 存在 严重 的 病态 特性 。 而 自前 通用 的 一 
些 LMI 求 解 工 具 ， 对 和 矩阵 的 条 件数 均 较 为 敏感 1， 故 这 类 病态 LMIs 不 适合 用 常规 LMI 工 具 
求解 。 对 稳定 性 分 析 问 题 如 此 ， 对 综合 问题 也 不 例外 。 
下 面 提出 时 澡 系 统 (2) 不 依赖 于 摄 动 参数 的 稳定 性 分 析 与 Hoo 滤波 器 设计 方法 。 


3 Ao MER BIRT 
假设 ， ARA (2) 是 渐进 稳定 的 。 
对 于 时 滞 系 统 (2)， 采 用 式 (4) 的 控制 规律 ， 设 计 五 - 滤波 器 。 
控制 器 规则 i: If £ (t) is Fii, and --- and €,(t) is Fy, then 
E(e)a(t) = Asi (t) + Briy(t), 
2(t) = Cyi&(t). 


(3) 
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利用 标准 模糊 推理 方法 ， 得 整体 控制 器 为 
E(e)è(t) = As (u)ê(t) + By(u)u(0), 


; (4) 
2(t) = Cs(u) a(t), 
其 中 
As(u) = So aids Br(p) = oe Cs(p = Smee 
t=1 i=l 
#1(t) 
E(t) = 
f | £2(t) 
21(t) € R”, ĉ2(t) CR 分别 为 滤波 器 的 慢 、 快 状态 向 量 。 
将 控制 器 (4) 应 用 到 系统 (2)， 得 到 闭环 系统 (5) 
Eei(e)&(t) = Aci(u)E(t) + Aga(p)E(t — 7(t)) + Ba(u)w(t), (0) = 0， (5) 


Z(t) = Cu(u)z(t), 
其 中 


B 
Aalu) = ane) : ， Balu) = 1H) | ; 
Bs(p)C2(p) Az (ue) Bs(p) Dai (ps) 


z(t) = | ce | ETE | b : | Egle) = | Ele) 0 | 


x(t — r(t)) 


&(t — r(t)) = | 
0 


EOS A= 20: 
定理 1 对 于 任意 给 定常 数 y > OMe e (0,s*] (0 < e* < 1)， 时 沾 模 糊 奇异 摄 动 系统 (2) 在 
在 一 个 y- 次 优 Hoo 滤波 器 (4) 当 且 仅 当 存在 矩阵 
可 Yu 0 
Yo. Y22 


X 0 
E 11 | 

X21 X22 
Xu E€ R”, Xa € R™*<™, Yu € R"Y",，Y22 € RT OPARE, IEEE 
MERE Qu e ROtmx(ntm), Ooo < ROtm)x(n+m) 和 适当 维 数 矩 阵 Mi, Ni, J， 使 得 矩阵 不 等 


A (6)-(9) 成 立 。 进 而 ， 如 果 和 矩阵 不 等 式 (6)-(9) 有 可 行 解 ， ae eee eee 
RH, 滤波 器 。( 下 面 式 中 * 号 表示 对 称 项 。) 


Ili < 0, i= 1,2, ,7, (6) 
My +0j <0, i<g<r, ij =1,2, 7 (7) 
Xıı — Yu > 0, (8) 


X22 一 Y22 > 0, (9) 


286 I E Kk 学 学 R 第 27 卷 


An * * * * * 
Aa A22 * * * * 
BEY BEX+ Di JF -I > * * 
aig Cii: — N; Cii 0 一 727 * * 
Adi Aii X 0 0 -G~A)Qu x 
0 0 0 0 0 -(1 — 8)Q22 
其 中 


ES d 


An = A Y +Y7A; + Qui + Ân, 
Ao = ATY + XTA;+ JiC2i + Mj +Q, 
A22 = ATX + XTA: + JjCzi + (JjCai)? + Qi, 


Agi=(¥-—X)'M;, Bri=(Y—X) I], Cr = Ni. (10) 
证 明 对 于 充分 小 的 e (0 <e <1), ELETA 
X: P 
P. 2 E 12 l 
Piz Le 
其 中 
X XZ P 0 Lı eLl, 
| Ap 2 aoa 
X21 X22 0 Piz La Lz 
AX € R"", Xo € R™™, Lu € R"™%", Ly € Rm 分别 为 对 称 正 定 矩 阵 ， X21 € 


R, La E R™™" AEREE, Pray E R°”, Prop E R™™ 为 对 称 和 矩阵 ， 则 存在 6 > 0, 
对 于 任意 的 es € (0, €5], HM P: WA Ba(e)P: = PI Eale) > 0. 


可 构造 Lyapunov 泛 函 


t 


V(&(t)) = 27 (t)Eci(e)Pe&(t) + f és ZT (0) QE(o)do > 0, 


-T 


其 中 


第 ?期 RA, AER: ARR Hoo R87 
Qu E ROM) Xintm), Qoy € Rintm)x(n+m) 分 别 为 对 称 正定 矩阵 。 
V(&(t)) = £7 (t)Ea(e)P-&(t) + £7 (t) Ea (€) P(t) 
+87 ()QE(t) — (1 — +E) (t — 7(t))QE(t — 7(t)) 
< T(t) [A] (u)Pe + PP Aalu) + QJ z(t) 
+87 (t—7(t))Ada(u)PeB(t) tw (t) BY (u)P.%(t) 
+27 (t)PP Aaa (u)E(t — T(t) + #7 (t)PF Ba(u)w(t) 
~(1 — B)ET(t — r(t))QE(t — r(t)) — w (t)w(t) + wT (t)w(t). 


CEASE RPI BM yT (t)2(t) (7 > 0 为 任意 给 定常 数 )， 然 后 对 于 任意 了 > 0， 不 等 式 
两 边 积分 得 


T 
V(#(T)) +7? f ZT (t)2(t)dt 


att) | [ ATOR + PZ Aalu) +Q4+77°CR(WCalu) o 
< f w(t) BT (u) Pe 一 7 * 
&(t — r(t)) Alet (Ut) Pe 0 -(1—£,)Q 
z(t) : 
w(t) dt + f wT (t)w(t)dt. 
0 
z(t — r(t)) 


V(a(T)) > 0。 因 此 ， 当 且 仅 当 书 满足 不 等 式 (11) 时 


T T Oo 
27 (t)z(t)dt < 7? f w? (t)w(t)dt < 7? Í wT (t)w(t)dt 
0 0 0 


成 立 。 我 们 有 
AT (u)P: + PP Aalu) +Q+Y ?CH(p)Ca(p) * * 
BT (p)P: -I * <0. (11) 
Agalu) Pe 0 —(1-8)Q 


显然 ， 不 等 式 (11) 可 改写 为 


AT (u)P + PT Aci(u) +Q+ 7 ?CH(u)Ca(u) * * 
BE (u)P 2 +0(e)<0, (12) 
Aal p)P 0 -(1-8)\Q 
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Xu ER?” Xy ER™™, Ly E€ 及 xm) Loz € 有 了 xm 为 对 称 正定 矩阵 ，Pi21 € R”*”, Pl22 € 
RM 为 对 称 和 矩阵 。 

根据 文献 [18,19] 中 的 结论 ， 可 得 存在 于 esi > 0， 对 于 任意 的 e e (0, ef], SAAS (13) RE 
时 ， 不 等 式 (12) 成 立 。 下 面 利用 


Al (u)P + PTAa(u) +Q +° (u)Calu) * * 
BT (u)P -I * <0, (13) 
AZalu)P 0 -(1—£8)Q 


对 不 等 式 (13) 应 用 Schur 补 定理 得 


Al(w)P+PTAi(H)+Q 4* * x 
BT 3 
cl (u)P I * * eS 0, (14) 
Calu) 0 -YI * 
Agalu) P 0 0 -(1-8)Q 
其 中 
oe | X Bp | 
Pio L 
令 
p= Z Sa 
Si2 S22 
则 
Z I I X 
P = l 
| Siz 0 | | 0 Pi2 | 
其 中 


Sin 0 
S12 = 
S122 S123 
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S121 € Rx”, S423 €eR™xm 为 任意 阵 。 
A 


Z I 
Ei = » = 
Si2 0 


| 
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对 不 等 式 (14) 左边 的 矩阵 分 别 左 乘 矩 阵 diag{ET, I, I, 1} MARIER diag{ E, I, I, 1}, 3 
结合 式 (5)， 得 


其 中 


邻 


其 中 


Way * * * 
Wo Woo * * 
BY (u) Woo -I + 


Ci(u)Z = Cilu) Si Cilu) 0 -yI 


Ad (1) 


0 


0 0 0 


Ai(u)X 0 0 ~(§-B)Qu 
0 —(1 ~ 6)Q22 


Way = ZAT (u) + A(u)Z + Z7QuZ + SHO22512, 


* 


* 


Wo) = AT (u) + XTA(U)Z + PBs (w)C2(u)Z + PRA) S2 + QuZ, 


Woo = A (u)X + X7A(u) + CF (u) BT (u) Pre + PEB CU) + Qui, 


Va2 


= BT (u)X + D5, (u) BF (u) Pi2. 


J(u) = PZB;(u), jp) = CHAS, 
SU FB BEAN AER (15) 可 改写 为 


Cu x fo. sp * * 
Pai Pes * oy x * 
BT (n) Pec S va * x 
C(u)Z- Ju) Cp) 0 -I * * 
Aq (H) AT(u)X 0 0 -(1-B8)Qi * 

0 0 0 0 0 —(1— B)Q22 


jn) = PhAs(u)Si2, 


Tuy = ZTA? (u) + Alu)Z + Z7 Qui Z + SEQ225h2, 


Tor = A? (u) + XT A(u)Z + I(u)Co(u)Z + J(u) + QuZ, 


T22 = AT (u)X + X7A(u) + CF (pw) JT (u) + J(u)C2(u) + Qu, 


T32 = By (u)X + Diy (u)I7 (u). 


(15) 


(16) 


(17) 
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对 不 等 式 (17) 左边 的 和 矩阵， 分 别 左 乘 矩阵 diag{2-T, I, I, I, I, 2-7 Si} MARA diag 
{Z7}, I, I, I, I, Si22-1}， 并 对 所 得 不 等 式 令 


Z =Y, M(u)=JSW)Y, Nw) =JW)Y, Qz =Y7SHQ2252Y, 


可 得 下 列 不 等 式 
On * * * * * 
日 ?1 O22 * * * * 
BE (HY oe e * <0, (18) 
Cilu)-N(u) Ce) 0 -yI * * 
43 (1) AT(u)X 0 0 -(1-8/)Qn * 
0 0 0 0 0 -(1 — B)Q22 
其 中 
ein = YTA(u) + AT (uY + Qui + Qo», 
O21 = AT(u)¥ + X7A(u) + J(u)C2(u) + M(H) + Qu, 
O22 = AT (p) X + XT A(u) + CZ HT (p) + J(u)C2(u) + Qu, 
932 = Bi (u)X + Da (u) J (p). 
由 式 (16) 得 


As(u) = P ÎS, Bile) = Pia J(u), Cpl) = J(u) Sia’, 
因此 ， 滤 波 器 (4) 的 传 函 矩阵 为 
Hz,(s) = JS [s1 — Pin” jp)Si2] Pia” J(u) 
= Ku)[s(I - XZ) - fp)] J) 


N(u)[sf — (Y — X) M (u) J Y - X) J(u), 


wa ‘ 
As(u) = (Y—X)"*M(n), Bri(p)= (YX) (uw), Cre) = N(u). 
结合 式 (2). (4), HASH (18) 改写 为 
> vil + 2. hiij (Iiz +My) <0, (19) 


t=1 
i<j 


因此 ， 存在 e > 0, e*= min{eġ, eï} 对 于 任意 的 e € (0,e*], 线性 矩阵 不 等 式 (6), (7) 成 立 。 
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另外 ， 由 Eale)Pe >0, E(e)X- = XT E(c) 可 得 
ET Ew(e)P EI 


ZTE(e)(XZ + Pi2S12)+ SBE(e)(Pi22Z+ LS12) (ZT XT + ST, Pi2)E(e) 
E(e)(XZ 十 Pl2S12) XTE(e) 


+O(e) > 0. (20) 


根据 PP-! = TAUCK [18,19] WE, MYA Cs > 0， 对 于 任意 的 e < (0e NRF 
式 (21) 成 立时 ， 不 等 式 (20) 成 立 。 注 意 到 


| ZTE(e) El) | 本 


(21) 
Ele) XTE(e) 


应 用 Schur 补 定理 ， 并 将 Y = ZRA: Yi 和 Yoo AM ARIE MM, HEER (8), (9) 成 
Ye 

因此 ， 存 在 s* > 0, e = min{eG, cf, 32}， 对 于 任意 的 e € (0,s*]， 闭 环 系统 (5) 满足 给 

E Ho 性 能 指标 7。 
at) | [ABWPHPIAAW+O i* z(t) 
z(t — r(t)) Ada (H) Pe —(1- BQ z(t — r(t)) 

根据 文献 [18,19] 的 结论 和 不 等 式 (14)， 得 V(z(t)) <0, AAA (5) HH BE. 

FEB ANSE SK (6)-(9) 是 与 e 无 关 的 线性 矩阵 不 等 式 组 ， 因 此 ， 可 以 应 用 Matlab 的 LMI 工 具 箱 
中 的 求解 器 feasp 求解 。 若 有 可 行 解 ， 根 据 定理 1 可 构造 系统 (2) 的 一 个 y- 次 优 五 滤波 器 。 考 
虑 求解 如 下 优化 问题 ， 如 果 该 问题 有 解 ， 则 结合 定理 1， 利 用 其 最 优 解 可 以 得 到 系统 (2) 的 最 
Th Hoo 滤波 器 。 本 文 应 用 LMI 中 的 求解 器 mincx 求解 。 


min y s.t. (6)-(9). (22) 


V(ž(t)) < | 


4 ”仿真 例子 


考虑 由 如 下 两 条 规则 描述 的 模糊 奇异 摄 动 系统 
R1: If zi(t) is Fi, then 


E(c)a(t) = Aiz(t) + Aaz(t — 7(t)) + Biiwt(t), 
z(t) = Cria(t), (23) 
y(t) = Caiz(b + Danw(t). 
R1: If zi(b is Fe, then 
E(e)a(t) = Aoa(t) + Aazz(t — 7(t)) + Bi2w(t), 
z(t) = C12z(t), (24) 


y(t) = C227(t) 十 D212w(t), 
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其 中 
一 0.1 10 一 4.6 10 一 0.1 0.2 
—] 一 1 一 ”一 1 一 0.1 一 0.1 
1 0 一 0.4 1.2 0 
so-| -| | y= a= l 
0 e —0.1 一 0.1 0.1 
Cu = C12 = C21 = C22 = [ 1 0 ] ， Dai = D212 =0.1, 7(t) = 0.5¢. 
对 应 模糊 集合 的 隶属 度 函数 为 


pi(zi(t)) = 0.5z1( 的 ， p2(z1(t)) = 1 — pi(z1(t)). (25) 


应 用 LMI 工 具 箱 及 其 求解 器 mincx， 求 解 式 (22)， 得 到 最 小 扰动 参数 Y = 4 时， 实例 系统 最 
Th Ho 滤波 器 参数 为 


1.0036 3.3433 _58.1831 -6.8660 0.5463 
Ap = ? 472 = 3 Br = 7 

2.2423 —25.6991 104.7622 —48.7340 1.1486 
9.0520 


Br = 
20.2301 


| , Cp = | -13364 -0.3322 ], Cpa = | -22.1426 -12.4140 | . 
以 初始 条 件 z(0) = [1,0,0,0], w(t) = sintrber-oat， 进 行 仿真 ， 图 1 至 图 4 为 对 应 < 
0.001, € = 0.1 时 的 闭环 系统 状态 响应 。 可 以 发 现 对 充分 小 的 <， 系 统 能 够 获得 满意 的 性 能 。 
注 : 图 1 至 图 4 中 的 xl,x2,xhatl 和 xhat2 分 别 表 示 闭 环 系统 状态 变量 zl (t), z2(t), ĉ1 (t), £2(t)。 


- €=0.001 时 x1 
£=0.1BYx1 


£70,001 BIx2 
5=0.1 时 Xx2 


t/s 


图 1: 闭环 系统 状态 变量 T1(t) 图 2: 闭环 系统 状态 变量 zz 人 (如 
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03 
=0.001 了 时 xhat1 E=0.001 时 xha 记 
5-0.1 时 xhat1 0.2 g-0.1 时 xhat2 
0.1 
0 
-0.1 
0 


图 3: ”闭环 系统 状态 变量 1(t) 图 4: 闭环 系统 状态 变量 2(t) 


结论 


本 文 建立 了 时 滞 模 糊 奇 异 摄 动 系统 模型 ， 通 过 Lyapunov 方 法 和 Schur 补 定理 ， 得 到 了 一 类 


时 滞 模 糊 奇 异 摄 动 系统 存在 模糊 五 。 滤波 器 的 充分 条 件 。 将 控制 器 设计 归结 为 与 摄 动 参数 无 
关 的 LMIs 求 解 ， 避 免 了 数值 求解 过 程 中 的 病态 问题 ， 并 能 够 获得 稳定 的 数值 解 。 对 充分 小 
的 =， 所 得 滤波 器 使 奢 环 时 滞 模 糊 奇异 摄 动 系统 渐进 稳定 ， 且 具有 给 定 的 Ho 性 能 指标 。 设 计 
方法 无 需 实 现 快慢 分 解 旦 适用 于 标准 和 非 标 准 非 线性 时 滞 奇 异 摄 动 系统 。 
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Hæ Filtering for Fuzzy Singularly Perturbed Systems with Time-delay 
CHEN Jin-xiang, YANG Wei-dong 


(College of Information Engineering, Beijing University of Science and Technology, Beijing 100083) 


Abstract: A modeling approach based on Takagi-Sugeno (T-S) is proposed and a fuzzy Hoo filter is 
developed for a class of nonlinear singularly perturbed systems with time-varying delay in states. The 
gains of the Ho filter are obtained by solving a set of the -independent linear matrix inequalities such 
that the ill-conditioned problem caused by the interaction of slow and fast dynamic modes is avoided 
effectively. The designed filter guarantees that the closed-loop fuzzy SPSs system is asymptotically 
stable and the prescribed level of disturbance attenuation is satisfied for sufficiently small e. The 
proposed approach can be applied to both standard and nonstandard nonlinear singularly perturbed 
systems. A simulation example is provided to illustrate the designed procedures. 

Keywords: singular perturbation; time delay; Ho filtering; linear matrix inequalities (LMIs); T-S 
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